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Abstract. It is shown that the nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) describes 
the shell-type nonthermal X-ray morphology of Cas A, obtained in Chandra observations, in a satisfactory way. The set of 
empirical parameters, like distance, source size and total energy release, is the same which reproduces the dynamical properties 
of the SNR and the spectral characteristics of the emission produced by CRs. The extremely small spatial scales in the observed 
morphological structures at hard X-ray energies are due to a large effective magnetic field 6 d ~ 500 //G in the interior which is 
at the same time not only required but also sufficient to achieve the excellent agreement between the spatially integrated radio 
and X-ray synchrotron spectra and their calculated form. The only reasonably thinkable condition for the production of such a 
large effective field strength is a very efficiently accelerated nuclear CR component. Therefore the Chandra data confirm first 
of all the inference that Cas A indeed accelerates nuclear CRs with the high efficiency required for Cas A to be considered as a 
member of the main class of Galactic CR sources and, secondly, that the nonlinear kinetic theory of CR acceleration in SNRs 
is a reliable method to determine the magnetic field value in SNRs. 
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1. Introduction 

Supernova remnants (SNRs) are the main sources of energy 
for the Interstellar Medium (ISM). They control the physical 
state of the ISM, and that presumably includes the nonther- 
mal component of Interstellar matter, often called the Galactic 
Cosmic Rays (CRs). Quantitative analysis of the SNR dynam- 
ics and its energetic particle production requires the determina- 
tion of the effective magnetic field in SNRs. For this we note 
that the magnetic field decisively influences the acceleration of 
CRs and their subsequent dynamics in SNRs. In the following 
we shall present evidence for the case of Cassiopeia A (Cas A) 
that the nonlinear kinetic theory for CR acceleration in SNRs 
(Berezhko et al. [T996l Berezhko & Volk [T997ll2000t does not 
only satisfactorily explain the observed SNR dynamics and the 
properties of the nonthermal emission, but that it also permits 
to determine the value of the effective magnetic field in SNRs, 
as already demonstrated for the case of SN 1006 (Berezhko et 
al. l2l)02ll2003al . 

The existence of a powerful population of relativistic par- 
ticles in SNRs is known from their nonthermal radiation. For 
some SNRs this emission is detected in a very wide range from 
the radio to the y-ray band. Cas A represents such an example. 



It is a shell-type SNR, and a bright source of synchrotron radi- 
ation observed from radio waves (e.g. Baars et al. !1977l > to hard 
X-rays (Allen et al. 119971 . These emissions are direct evidence 
for the existence of a large number of relativistic electrons pre- 
sumably accelerated at the shock, with energies reaching about 
10 TeV. Direct experimental knowledge of nuclear CR produc- 
tion can only be obtained by high energy y-ray measurements. 
If protons are efficiently accelerated in Cas A, then the n°- 
decay y-ray spectrum as a result of inelastic collisions with 
the background nuclei must extend beyond 1 TeV with a hard 
power-law spectrum. The detection of a corresponding if weak 
signal in TeV y-rays has been indeed reported by the HEGRA 
collaboration (Aharonian et al. 1200 11 1. 



Send offprint requests to: H. J. Volk 



According to the analysis by Borkowski et al. ( 1996) which 
is consistent with the observed dynamics of Cas A, the super- 
nova shock expands into a non-uniform circumstellar medium 
strongly modified by the intense progenitor star wind. This cir- 
cumstellar environment consists of a tenuous inner bubble, cre- 
ated by the Wolf-Rayet wind, a dense shell of swept-up slow 
Red Supergiant (RSG) wind material, and a subsequent unper- 
turbed RSG wind. The medium further out is not relevant here. 
Since the shock has akeady passed through the dense shell 
and has swept up a mass of about 5M (Favata et al. 1997 1, 
a detectable 7r°-decay y-ray flux is to be expected. Recently 
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Chevalier & Oishi (2003 1 have argued that a Wolf-Rayet phase 
and the corresponding modification of the preexisting RSG 
wind is not required. To reach a final conclusion about this 
important issue a detailed dynamical and kinetic consideration 
has to be made, similar to the one contained in our previous 
paper (Berezhko et al. 2003b I. At the moment we can only say 
that the global characteristics of nonthermal radiation of Cas A 
are expected to be not sensitive to these details of the struc- 
ture of circumstellar medium, because they are mainly deter- 
mined by the total amount of the swept up gas and the mag- 
netic field value in the remnant. Since the observed hard X-ray 
synchrotron emission of Cas A comes from the thin region near 
the current SN shock position, the properties of this radiation 
are not expected to be very sensitive to the previous history of 
the SNR evolution. Here we use our model (Berezhko et al. 
2003b I which agrees with all previous data, leaving the ques- 
tion about the role of Wolf-Rayet stage for future scrutiny, and 
test its consistency with the hard X-ray morphology. 

Comparing the nonthermal radio and X-ray synchrotron 
data with the calculated spectrum of the energetic electrons, 
Berezhko et al. ( 2003b I have inferred that the existing data re- 
quire very efficient acceleration of CR nuclei at the SNR blast 
wave which has converted already several percent of the initial 
SNR energy content into nuclear CR energy despite the neces- 
sity of renormalization of the overall nuclear CR content due 
to the nonuniform injection of suprathermal ions (Volk et al. 
2003 1. At the same time a large interior magnetic field strength 
Z?d ~ 500//G is required. The large value of B d is expected and 
possible in this case, using the spatially integrated radio and X- 
ray synchrotron spectra to empirically determine the effective 
magnetic field strength and the injection rate in the theory (see 
also Volk l2003l . 

In fact there are very recent measurements of details of the 
X-ray morphology of Cas A with the Chandra telescope by 
Vink & Laming (2003). These authors have used the detailed 
Chandra morphology to estimate the magnetic field strength 
by assuming the small spatial X-ray scale to be determined by 
postshock synchrotron losses. Our aim is to show that the ac- 
tual B -field strength coincides with that inferred from the radio 
spectrum and nonlinear acceleration theory. We confirm that 
the field is indeed amplified in Cas A, although it turns out that 
the actual value is significantly larger than estimated by Vink 
& Laming (E003- 

2. Results and discussion 

Compared with SN 1006 which exploded into an essentially 
uniform ISM, Cas A is a very different and much more complex 
object. The dynamics of the SN shock and the major fraction of 
the emission produced by accelerated CRs are predominantly 
determined by the parameters of the swept-up RSG shell. They 
depend less sensitively on the parameters of the medium near 
the current shock position since that presumably consists of 
free RSG wind material and contains little mass. Nevertheless, 
since the required magnetic field in the interior downstream of 
the shock is so large that it strongly influences the energy spec- 
trum of the CR electrons as a result of their synchrotron losses 
(Berezhko et al. 2003b ), the hard X-ray synchrotron emission 



should come predominantly from the relatively thin postshock 
region, because deeper inside the remnant the upper cutoff en- 
ergy of CR electrons becomes too low. Therefore one would 
expect that the nonlinear kinetic theory should not only repro- 
duce the spectral characteristics of Cas A in a satisfactory way, 
but should also be able to explain the observed fine structure of 
the nonthermal emission in hard X-rays. 

For the comparison with the Chandra data we shall not 
present the details of our model. They have already been de- 
scribed in the above paper (Berezhko et al. 2003b I. In order 
to give a clear interpretation and simple estimate of the spa- 
tial variation of the synchrotron emission near the SNR shock 
we shall use some simple analytical approximations, extending 
the discussion of SN 1006 (Berezhko et al. 2003a}. As a second 
step we shall make a comparison of the Chandra data with the 
results of the nonlinear theory. 

It was shown (Berezhko et al. 2003a) that due to the strong 
synchrotron losses the very energetic electrons occupy only a 
thin region behind the shock of thickness 

l 2 X = (2/ d2 )- 1 ( Vl+4/d2/^2 - 1), (1) 

where Uip) = k(p)/u and l c (p) - r(p)u are the diffusive and 
convective length scales respectively, r(p) = 9m 2 c 2 /(4r^B 2 p), 
is the synchrotron loss time, m e is the electron mass, ro denotes 
the classical electron radius, B is the effective magnetic field 
strength, <e(p) = psv/3 is CR diffusion coefficient, and Pb(p), 
p and v are the particle gyroradius, momentum and velocity, 
respectively. In addition 112 = V s /cr is the downstream plasma 
speed relative to the shock front, V s is the shock speed, and cr 
is the total shock compression ratio. 

Depending upon the ratio /d/Z c we have two different ex- 
treme cases. For electron momenta p corresponding to l c » 41& 
(weak losses) we have I2 = k- In this case the loss time r is 
much larger than the acceleration time r acc ~ k/u 2 and there- 
fore the losses do not influence the electron spectrum during 
their acceleration. The electron spectrum has a power law char- 
acter up to a cutoff momentum p max > p, which exceeds the 
considered momentum p. Later on, the losses in the down- 
stream region lead to a decrease of the cutoff momentum so 
that at the distance r — R s - l c it drops to the value p max (r) = p. 
Here R s is the shock radius. In the opposite extreme case of 
strong losses l c <K 4/d, we have I2 = y[i<2J2- which is in fact 
independent of momentum if particle diffusion proceeds in the 
Bohm limit. In this case losses become significant already dur- 
ing the electron acceleration and therefore the cutoff momen- 
tum p m . dx is lower than p, so that the corresponding electrons 
belong to the exponential tail of the distribution. Since the non- 
thermal X-ray emission of Cas A belongs to the steep tail of 
the synchrotron spectrum, Eq.Q suggests that / c <K Zd for the 
radiating electrons. 

Since the synchrotron data are given in terms of the emis- 
sion frequency v, it is useful to rewrite Eq.([0, taking into ac- 
count the approximate relation v oc Bp 2 : 

h = yfwi/i Vl + 6 2 - 6), (2) 
where 

5 1 = l C 2/(4l d2 ) = 0A2[c/(r v)][VJ(crc)] 2 . (3) 
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We then have the situation that the synchrotron emissivity 
q v (e v , r) is everywhere close to zero except for a thin radial re- 
gion of thickness Ar ~ 1% just behind the SN shock, since the 
emissivity from the upstream region r > R s is insignificant. 

In projection along the line of sight, the radial emissivity 
profile determines the remnant's surface brightness J v . For the 
X-ray energy e v it has the form 



J v (e v ,p)-2 J dxq v (e v ,r- ^Jp 2 

J a 



X 2 ), 



(4) 



where p is the distance between the center of the remnant 
and the line of sight. In the case of a spherical shock a = 

To estimate the shape of the expected profile J v (p) one can 
consider the simple case, when the emissivity has a form 



q v (e v , r) = q 2 (e v ) exp[(r - R s )/l 2 ] 



(5) 



for r < R s . In the case of interest, I2 <s /? s , the main contri- 
bution to the integral @ comes from the thin regions near the 
boundary x = a. To perform the integration along the line of 
sight one can therefore expand the quantity r = yp 2 + x 2 in 
Eq.l|5j in powers of x - a up to the second order term which 
results in: 



2q 2 Rsk 



■p- 



R 2 * - V 



1 



R 2 s -2p 2 (z + 2)/2p 2 -R 



R 2 z 



r-: 



(6) 



where z = (R 2 - p 2 )/(R s l2)- This approximate expression is ac- 
curate to about one percent for the case of a small thickness of 
the emission layer l 2 <k Q.IR S . As it is clear from this expres- 
sion the distant observer will see the profile J v (p) which, going 
from large distances to the center of the SNR, starts at p = R s , 

reaches its peak value 7 max = 2q 2 Jr 2 — p^ at p m ~ R s - 0.85/2 
and smoothly drops to / m i n = 2/2^2 at the SNR center p = 0. 
Therefore the outer scale L\ ss I2/2 characterizes the emission 
profile at p > p m , whereas for the inner scale, which charac- 
terizes the profile at p < p m , we have L2 ~ 13Li. This means 
that the width L — L\ + L2 ~ lh of the observed brightness 
profile J v (p) is always appreciably wider than the width l 2 of 
the emissivity layer q v (r), simply due to the geometry of the 
system. 

According to @ the downstream magnetic field is 



B d = [3m 2 e c 4 /(4er 2 l 2 )] l/ \ Vl +( 5 2 - g) 



-2/3 



(7) 



The case of strong losses l c <K 4Zd corresponds to 6 = 0. 
According to Vink & Laming (2003) the total thickness of 
the observed profile at e v > 2.7 keV is as small as 1.5". 
However, following the above consideration, only the small 
fraction 1.5"/7 w 0.22" corresponds to the thickness of the 
downstream emission region h - For a distance d — 3.4 kpc and 
an angular radius of 2.6' for Cas A this gives h ~ 1 . 1 x 10 16 cm. 
In the case of strong losses (6 = 0) B& « 470 ^G, independently 
of synchrotron energy and shock speed. For Vs = 3000 km/s 
and cr = 6, corresponding to the overall nonlinear model 



1 — 1 — 1 — 1 — 1 — 1 — r_ 
e, >2.7 keV 




ip (arcsec) 

Fig. 1. The X-ray emissivity for X-ray energies 2.7 < e v < 
9 keV as a function of angular radial distance. The Chandra 
measurements are shown (Vink & Laming 2003 1. 



(Berezhko et al. 2003b 1, 5 2 = 0.054 resulting in B d * 550 pG. 
For a larger shock speed V s = 5000 km/s, as derived from X-ray 
measurements and used by Vink & Laming, we have 6 2 = 0.15 
and B d ~ 610 pG. 

The corresponding magnetic field value is by a factor of 



about five larger than estimated by Vink & Laming (2003 1. 
There are two reasons for such a large discrepancy (i) only a 
small part of the observed profile width L corresponds to the 
thickness I2 of the emission region (ii) the correct formula Q, 
cf. Volk et al. (1981), implies l C 2 < h- Using l C 2 = L, as was 
done by Vink & Laming, therefore doubly underestimates B&, 
since l C 2 °c B^ 2 . 

The numerically calculated brightness profile J v (p) = 

r-OS 

J de v J v {e v , p) for the X-ray energy e v > e\ — 2.7 keV is shown 
in FigtO The sharpest experimental X-ray brightness profiles 
obtained by the Chandra observers (Vink & Laming 2003 1 are 
shown as a function of angular distance i/r = 2.6x60(p-R s )/Rs 
in arcseconds, taking into account that the angular radial size 
of Cas A is 2.6'. Since the absolute values of the measure- 
ments are not known, the theoretical profile is normalized so 
as to obtain the best fit to the experimental profile. The very 
sharp part on the right side of the calculated profile J v {\jj) cor- 
responds to the shock front position p - R s . Compared with 
Vink & Laming ( 2003 1 all experimental points are shifted by 
Ai// = 0.3" to the left. One can see that the experimental values 
agree quite well with our calculations, which yielded a post- 
shock value of the magnetic field strength at the current evolu- 
tionary epoch (see Berezhko et al. 2003b 1 of Z?d = 480 pG. 
Taking into account that the Chandra angular resolution is 
0.5" we conclude that the actual downstream magnetic field 
in Cas A is indeed Bd ~ 500 pG with an uncertainty of about 
30% which encompasses the above simple estimates of the lo- 
cal field. 

There are a number of reasons for a broadening of the ob- 
served profile, compared with the ideal case of a spherical shell. 
First of all, it is clear from the observations that the actual shock 
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front deviates from a strictly spherical form. The wavy shape 
of the shock front can be produced as the result of small scale 
density inhomogeneities of the ambient ISM. Any small-scale 
distortion of the spherical emission shell leads to a broaden- 
ing of the observed brightness profile. This explains why the 
observed X-ray radial profile is slightly broader than the the- 
oretical profile (see Fig.l). Therefore it is hard to believe that 
the actual magnetic field in Cas A could be significantly lower 
than the above value. 

Since at the current epoch the supernova shock propa- 
gates through the free RSG wind, where the required effec- 
tive magnetic field would be so strong that the Alfven speed 
is c a m 160 km/s (Berezhko et al. 2003b i, this field exceeds 
the preexisting wind field by at least an order of magnitude, 
because the wind should presumably be superalfvenic, that is 
c a < V w , where V w ~ 10 km/s is the RSG wind velocity. From 
this argument, the required strong downstream magnetic field is 
inevitably the result of considerable amplification near the su- 
pernova shock. As argued for the case of SN 1006 (Berezhko 
et al. I2003a> such a strong magnetic field amplification can 
only be produced nonlinearly by a very efficiently accelerated 
nuclear CR component. In this case their number, consistent 
with all existing data, is so high that they are able to strongly 
excite magnetohydrodynamic waves, and thus to amplify the 
magnetic field, and at the same time to permit efficient CR 
scattering on all scales, approaching the Bohm limit (Bell & 
Lucek 2001). The same large effective magnetic field is re- 
quired by the comparison of our selfconsistent theory with the 
synchrotron observations. 

For such a high downstream magnetic field * 0.5 mG, 
which leads to strong synchrotron losses of the CR elec- 
tron component, the leptonic contribution to the TeV y-ray 
emission is so small that it is far below the detected HEGRA 
flux (Berezhko et al. l2003bt . 
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3. Summary 

We conclude that the Chandra data confirm the magnetic field 
amplification and the very efficient acceleration of nuclear CRs 
in Cas A, predicted by Berezhko et al. ( 2003b I. This efficiency 
is consistent with the requirements for the Galactic CR energy 
budget. At the same time, together with SN 1006 (Berezhko et 
al. 2002; 2003a I this is now a second - and astrophysically very 
different - case, where the prediction of the nonlinear kinetic 
model for CR acceleration in SNRs is fully confirmed exper- 
imentally. Therefore we conclude that this theory is not only 
able to describe the CR dynamics and acceleration in SNRs, 
but that it constitutes in addition a reliable method to quantita- 
tively determine the effective magnetic field strength which is 
produced in the acceleration process. 

Since our model uses the Bohm diffusion of CRs, which is 
consistent with the field amplification picture, one can consider 
the consistency of the magnetic field values required to explain 
the global properties of the synchrotron emission on the one 
hand (Berezhko et al. 2003b I, and to reproduce the fine struc- 
ture of the X-ray radiation on the other as observational support 
for the Bohm limit of CR diffusion near the strong SNR shock. 



